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We propose that the observed spin-relaxation in bilayer graphene is due to resonant scattering by 
magnetic impurities. We analyze a resonant scattering model due to adatoms on both dimer and 
non-dimer sites, finding that only the former give narrow resonances at the charge neutrality point. 

Opposite to single-layer graphene, the measured spin-relaxation rate in graphene bilayer increases 
with carrier density. Although it has been commonly argued that a different mechanism must be at 
play for the two structures, our model explains this behavior rather naturally in terms of different 
broadening scales for the same underlying resonant processes. Not only our results—using robust 
and first-principles inspired parameters—agree with experiment, they also predict an experimentally 
testable sharp decrease of the spin-relaxation rate at high carrier densities. 
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Understanding spin-relaxation is essential for de¬ 
signing spintronics devices nua. Unfortunately, spin- 
relaxation in graphene structures has been a baffling 
problem [3]. While experiments in both single layer 
graphene (SLG) [414TT] and bilayer graphene (BLG) [7118] 
yield spin lifetimes on the 100 — 1000 ps time scale 
(the highest values achieved in graphene/h-BN struc¬ 
tures Hang), theories based on realistic spin-orbit cou¬ 
pling and transport parameters predict lifetimes on the 
order of microseconds mill. 

While the magnitudes of the spin-relaxation rates of 
SLG and BLG are similar, the dependence of the rates 
on the electron density is opposite in the two systems. In 
SLG the spin-relaxation rate decreases with increasing 
the carrier density EHS], in BLG the spin-relaxation rate 
increases [aiH]. Since the diffusivity in the investigated 
samples decreases with increasing the electron density, 
it has been a common practice to assign two different 
mechanisms to both structures: the Elliott-Yafet mech¬ 
anism [25l [26] to SLG O H |9l [10] and Dyakonov-Perel 
mechanism [27] to BLG [7]-[9]. 

The main problem with that assignment is quantita¬ 
tive. Spin-orbit coupling in graphene [28] is too weak to 
yield such a small spin-relaxation time. An explicit first- 
principles calculation [23] predicts that one would need 
0.1 % of adatoms to give 100 ps spin lifetime. Recently 
a new mechanism for SLG was proposed m (see also 
Ref. [30]), based on resonant scattering off local mag¬ 
netic moments. It gives the observed spin-relaxation 
times with as little as 1 ppm of local magnetic mo¬ 
ments and also agrees with the experimental behavior 
for SLG of decreasing the spin-relaxation rate with in¬ 
creasing electron density. Where do these local moments 
come from? It was theoretically predicted that adatoms 
such as hydrogen ISBEg, but also chemisorbed organic 
molecules [33] can be responsible. Experimentally it was 


demonstrated that hydrogen adatoms indeed induce lo¬ 
cal moments [34[ [35] , but even untreated graphene flakes 
were shown to exhibit 20 ppm spin 1/2 paramagnetic 
moments [36]. The most natural candidates for reso¬ 
nant magnetic scatterers appear to be polymer residues 
from different fabrication steps of graphene devices. Our 
mechanism is also in line with mesoscopic transport ex¬ 
periments naEi] which found a strong evidence for the 
local magnetic moments in the dephasing rates. Our the¬ 
ory does not work for high adatom concentrations (say, 
above 0.1%), at which the induced magnetic moments 
seem to form a fluctuating magnetic-field network [34] . 

In this letter we propose that the spin-relaxation in 
BLG is caused by the same mechanism of resonant mag¬ 
netic scatterers. We show that (i) adatoms (we model 
specifically hydrogen) on dimer and non-dimer sites of 
BLG give different resonance energies and resonance 
widths, (ii) the calculated spin-relaxation times are in 
quantitative agreement with experiment, (iii) the oppo¬ 
site trends of the spin-relaxation rate in SLG and BLG 
are due to different scales of the energy fluctuations 
(caused by electron-hole puddles) in the two structures, 
reflecting their different density of states (DOS), (iv) the 
spin-relaxation rate in BLG should reverse its trend and 
decrease with increasing electron density at high densi¬ 
ties, making an experimentally verifiable prediction. As 
in SLG, also in BLG resonant magnetic scatterers are 
spin hot spots [39]: affecting spin but not momentum 
relaxation. 


Model Hamiltonian. We consider a single adatom on 
AB stacked bilayer graphene sitting on either a dimer 
or a non-dimer position. The full model Hamiltonian is 
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FIG. 1. (Color online) Calculated electronic structure of bilayer graphene with hydrogen adatoms. Panels (a) and (b) are 
for dimer adatoms, (c) and (d) for non-dimer ones. In (a) and (c) we plot the electronic band structures: dotted lines are 
spin-unpolarized first-principles calculations using a 7x7 supercell, while solid lines are tight-binding fits as described in text. 
Panels (b) and (d) show unperturbed, Qq (E) + Qq (E), and perturbed, lZc{E) (with adatom concentration of 77 = 0.05%), DOS 
per atom and spin. Dimer adatoms (b) show a narrow resonant peak near the charge neutrality point at E^es — 22.5 meV 
with the full width at half maximum P ~ 8.4 meV. Non-dimer adatoms (d) induce a broad resonance at E^es — 26.1 meV 
with r ~ 165.2 meV. For plotting DOS we perform running averages of 20 meV. Insets: schemes of the tight-binding model 
Hamiltonian, i7o + 77', Eqs. Q and 


Hq -f 77', where 

Ho = -to \a^^) {b^^\+tiY\(La) {bma\+^-C; 

{m,n)a ma 

A e {t ,b} 

( 1 ) 

is the unperturbed BLG Hamiltonian with intralayer 
nearest-neighbor hopping to = 2.6 eV, and direct inter¬ 
layer hopping ti = 0.34 eV [40]. We neglect indirect in¬ 
terlayer parameters which yield fine features of the energy 
bands (such as warping and electron-hole asymmetry) as 
unimportant for our purposes. The first sum runs over 
(m, n) nearest neighbors in the top (A = t) and bottom 
(A = b) layers. The second sum runs over the a-sites of 
the A sublattice of the top layer and 6-sites of the B sub¬ 
lattice of the bottom layer. State |c^cr) represents the 
spin a carbon 2pz orbital on sublattice c = {a, 6} and 
site m in layer A. The eigenstates of 77o are miiig 

£^J(k) = I (Mb + , (2) 

where the index a labels the conduction {a = +) and 
valence {a = —) bands, and /i stands for the high {p = +) 
and low (/i = —) energy bands with respect to the charge 
neutrality point, for /(k) see [43]. 

We place the adatom on the top layer either on a dimer 
site Cd, which is an site, or on a nondimer site Cnd^ 
which is a 6^ site. The adatom also carries a local mag¬ 
netic moment coupled with the electron spins via the 
exchange interaction J. The Hamiltonian 77' is [29] . 

H’ = Y^\ha) {K\+uj{\K) {c,\ + \c^) (3) 


where \hcj) is the adatom orbital with spin a. This or¬ 
bital has on-site energy 5 and is connected to the site C 
on the bilayer with hopping energy uj. The spin operators 
s and S, which are the Pauli matrices in the correspond¬ 
ing spinor spaces, are for itinerant electrons and local 
magnetic moments, respectively. 

To obtain realistic parameters for the adatom Hamilto¬ 
nian 77', we performed first-principles calculations with 
Quantum ESPRESSO [44] using a 7 x 7 graphene su¬ 
percell with a single hydrogen adatom. In agreement 
with previous studies [45] we found that hydrogen on 
both dimer and non-dimer sites induces local magnetic 
moments of 1 Bohr magneton per unit cell. However, for 
fitting the orbital parameters of 77', namely e and cj, we 
used the spin-unpolarized first-principles band structure 
and set J = 0 in the tight-binding calculation. For the 
dimer site we select 5 = 0.25 eV and uj = 6.5 eV, while 
for the non-dimer one e = 0.35 eV and uj = 5.5 eV. Fig¬ 
ures Qa) and (c) show that the fits are very good. How¬ 
ever, these fitted parameters are not unique, as a larger 
neighborhood of values offers a comfortable agreement 
with first-principles data. We have selected the values 
which are close to the uniquely fitted SLG orbital pa¬ 
rameters. For the exchange coupling we take the same 
(unbiased) value as for SLG [29], J = —0.4eV. 

Resonant scattering. We transform the adatom 
Hamiltonian 77' into the singlet (7 = 0) and triplet (7 = 
1) basis \ci^rni) = |c) G) |7, m^) [label m£ = —7,..., 7 is the 
total spin projection] and eliminate the adatom orbital 
\h) by downfolding. This gives the energy dependent 
perturbation H'{E) = 'Z,e,me allow¬ 

ing us to analytically calculate the T-matrix, T(£') = 
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FIG. 2. (Color online) Calculated spin-relaxation rates 1 /ts as a function of energy (carrier-density) for impurity concentration 
rj = O.lTppm. (a) Zero temperature, unbroadened, rates for dimer (red solid) and non-dimer (black solid) adatoms, as well as 
the resulting average 1/ts (blue dotted). For reference the SLG calculation is also shown (gray dashed-dotted), (b)-(d) Spin- 
relaxation rates for three representative temperatures. Theoretical data (blue solid) are broadened, simulating the presence 
of electron-hole puddles, with a Fermi level smearing of 23meV. Circles and diamonds represent data points from Aachen- 
Singapore (AS) [7] and Riverside (R) [8] experiments, respectively. The two shoulders (spin-relaxation edges) at ±100 meV are 
exchange-split resonances. At high carrier densities the model predicts a decrease of the spin-relaxation rate. 


\ce,me) (Q.mJ, where 

Ve{E) = = i-Ve(E)Gc{E) ■ 

Here Gc{E) = Ac{E) — muc^E) is the C-site projected 
Green’s function per atom and spin of the unperturbed 
BLG with 



At=± 

where D = \/ y/Sirto — 6eV is the effective BLG band¬ 
width and Ac equals zero for Cd-site and one for Cnd-site, 
respectively. 

We first analyze orbital resonances of H' (set J = 
0) by plotting in Fig. Bb) and (d) the perturbed 
DOS per atom and spin, lZc{E) = X]/x=±^o(^) “ 
(77/7r) Im { [— ^Gc(L^)] T^(±^, J = 0 )}, where 77 is the 
adatom concentration per carbon atom and Qq{E) = 
(2|F;| - /iti)/(4D2) 0(D - \E\) 0(|F;| - /iti) is the un¬ 
perturbed bilayer DOS per atom and spin for the high 
(/i = ±) and low (// = —) energy band, for details 
see [43]. As seen from Fig. [^b), the dimer site hydrogen 
chemisorption induces a pronounced narrow resonance 
near the charge neutrality point. In contrast, non-dimer 
adatoms, see Fig. [^d), give a broad and shallow reso¬ 
nance. This striking difference is explained by consider¬ 
ing the character of the resonance states. In a monolayer 


graphene an adatom on the A site induces a resonance 
state which is localized mainly on B sublattice. Thus, an 
adatom on a dimer site induces a resonance state which 
is spread mainly on the non-dimer sublattice and hy¬ 
bridizes only little with the other layer, keeping the res¬ 
onance narrow. If the adatom is on a non-dimer site, the 
resonance state is mainly on the dimer sublattice which 
couples to the other layer, causing a leakage of the state 
and broadening of the resonance. The same behavior is 
seen in vacancy models [46]. 

Spin-flip scattering and spin-relaxation rate. The T- 
matrix allows us to compute the spin-flip rate for a sin¬ 
gle scattering event by adatom at site C (dimer or non¬ 
dimer) |k^(F^),t) ^ \q^{E'),l) between bands fi and 

1/ im, 

mEWUE') S{E - E') , (7) 

where the site and band dependent projections Fq{E) = 
2{\E\-^xAch)/{2\E\-^it^) e{D-\E\) e{\E\-fih), see 
also m- The exchange-induced spin-flip function is 


/b(^) 


^ {£-l)Ve{E) 

,^p-Vi{E)GciE) 


(8) 


The spin-flip rate does not depend on the relative ori¬ 
entation of k and q, since the energy dispersion in our 
model has rotational symmetry. However, the spin-flip 
rate is very different for dimer and non-dimer adatoms. 
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To obtain the spin-relaxation rate l/r^ we sum over 
different partial rates and obtain 

1 2»,c ,p,K(-EW(fi) + P5(£)eo(£)]" 

where the labels + and — denote BLG high and low en¬ 
ergy bands entering the definitions of Pq and Qq given 
in the text. To get the final spin-relaxation rate we take 
an unbiased average over the dimer and non-dimer sites, 
1/ts = 1 / 1 ). This is plotted in Fig. |^a) 

and compared with SLG. Two pronounced shoulders—we 
call them spin-relaxation edges —in emerge from 

the exchange splitting of the orbital resonance seen in 
DOS at Fig. [^b), just like for SLG, although the peaks 
in BLG are more separated due to the energy renormal¬ 
ization by the interlayer coupling. In contrast, non-dimer 
adatoms show a rather flat behavior with respect to the 
energy, reflecting the broad resonance of the perturbed 
DOS, Fig. [^d). Non-dimer adatoms still induce a large 
I/ts since they strongly perturb the low-energy states 
which are localized on the non-dimer sites. This behav¬ 
ior is encoded in the low-energy site projection P^{E), 
in Eq. which is at low energies much larger for non¬ 
dimer than for dimer adatoms. 

Comparison with experiments and eontrasting single 
and hilayer graphene. Gomparison with experiments re¬ 
quires temperature and electron-hole puddles broadening 
of I/ts. Temperature broadening is due to population 
smearing, {—dfo/dE), where /o is the Fermi-Dirac dis¬ 
tribution. The puddle broadening is modeled as a con¬ 
volution with a Gaussian kernel of standard deviation 
cTbr- For bilayer we use abr — 23meV, which corre¬ 
sponds to density fluctuations An of 8.5 x 10^^/cm^, 
following experimental estimates [48]. In Figs, [^b)- 
(d) we present the main results of this paper, the fully 
broadened spin-relaxation rates compared with Aachen- 
Singapore (AS) [7j and Riverside (R) |8| experiments. 
Glearly the two experiments are somewhat at odds, but 
they display consistent behavior at low temperatures. We 
adjust the local moment concentration to describe the 
AS data, obtaining r] = 0.17ppm. All other parameters 
are as obtained from the orbital fits. The agreement at 
low temperatures is especially good. At high tempera¬ 
tures the overall shapes differ, but the two experiments 
differ as well. This experimental discrepancy further un¬ 
derlines the extrinsic character of the spin-relaxation in 
BLG. It is likely that the relative population of dimer 
and non-dimer adatoms changes with temperature, dif¬ 
ferently in different samples, reflecting the idiosyncrasy 
of the experimental data. However, our calculation gives 
a rather robust prediction at low temperatures: at high 
earner densities, above the spin-relaxation edge at about 
5 X 10^^/cm^, the spin-relaxation rates should start to 
deerease. 

Perhaps the most pressing remaining question is: 
Given the same resonant spin-relaxation mechanism for 



FIG. 3. (Color online) The effect of electron-hole puddles 
on spin-relaxation in SLG and BLG. From top to bottom: 
The spin-relaxation rate exhibits two resonance peaks due to 
singlet-triplet splitting AFst- The splitting of the peaks is 
greater in BLG. The peaks are broadened by temperature and 
carrier density fluctuations An which is very different for SLG 
and BLG, due to their different DOS. For a given temperature 
and density fluctuation An the energy smearing in SLG crbr — 
AFst, while in BLG Cbr ^ AEst- After broadening the spin- 
relaxation rate around the charge neutrality point in SLG 
has the opposite trend as the unbroadened rate. In BLG the 
original trend is preserved. 


single and bilayer graphene, why do their spin-relaxation 
rates have the opposite trends as functions of charge den¬ 
sity mi? Our mechanism offers a natural, and perhaps 
mundane answer: electron-hole puddles. At low tem¬ 
peratures and in the absence of density fluctuations the 
two structures should exhibit the same trend, namely, 
an increase of the spin-relaxation rate going away from 
the charge neutrality point. In SLG the behavior is ex¬ 
actly opposite. The reason is offered in Fig. In SLG 
the carrier density fluctuations lead to a large Fermi 
energy smearing (dbr = 91 meV versus dbr = 23 meV 
in BLG for the same carrier density fluctuation An of 
8.5 X 10^^/cm^). Averaging over the Fermi energy of the 
singlet-triplet split spin-relaxation peaks then inverts the 
shape of the spin-relaxation rate around the Dirac point. 
In bilayer, due to its greater density of states, the energy 
broadening is much more modest, and the experiments 
(unless their samples would exhibit large variations of the 
electronic densities) find the behavior as expected for an 
unbroadened system. Figure also shows the origin of 
the spin-relaxation edge and the robustness of our predic¬ 
tion of the decrease of the spin-relaxation rate at greater 
electron densities. At high temperatures (above 100 K), 
it is enough to invoke thermal broadening to see the trend 
reversal even in ultraclean SLG, as its resonance peaks 
are closer than those in BLG m- The picture given in 
Fig-i could be used to analyze experimental trends in 
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spin-relaxation in both SLG and BLG. 

In conclusion, we showed that resonant scattering by 
local magnetic moments quantitatively accounts for the 
experimental data. This spin-relaxation mechanism also 
explains the apparently striking opposite behavior of the 
measured spin-relaxation of SLG and BLG, offering a real 
alternative to quantitatively unsubstantiated but often 
made assignment of the two distinct trends as Elliott- 
Yafet and Dyakonov-Perch Finally, our model makes a 
specific prediction of reversing the increase of the spin- 
relaxation rate in graphene bilayer with increasing carrier 
density, at high densities, accessible experimentally. 

We thank B. Beschoten and R. Kawakami for provid¬ 
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ported by DFG SFB 689 and GRK 1570, and by the 
FU Seventh Framework Programme under Grant Agree¬ 
ment No. 604391 Graphene Flagship. 
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FIRST-PRINCIPLES METHODS 

We performed first-principles calculations with Quan¬ 
tum ESPRESSO [1] for hydrogen adatoms on BLG, 
using 7x7 supercells, in the slab geometry with vac¬ 
uum spacing of 15 A. The atomic positions were re¬ 
laxed using the quasi-newton algorithm based on the 
trust radius procedure. The exchange-correlation poten¬ 
tial was treated within the generalized gradient approx¬ 
imation [2], and the van der Waals interaction was in¬ 
cluded within a semiempirical approach [3, 4]. We used 
nltrasoft psendopotentials [5] with plane wave and charge 
cutoff of 30 Ry and 300 Ry, respectively. 

ELECTRONIC STATES OF BILAYER 
GRAPHENE 

Direct diagonalization of Hamiltonian Hq of the 
manuscript in the associated Bloch basis, 

'' m 


[E+-si^(k)] \ AsUenters 
the definition of the T-matrix, T = E/(l — GcV) |c)(c|, 
it specifies resonant energies of the system at which 
Gc(Eres) 1/E. In the long-wavelength limit with en¬ 
ergy cut-off D, the Green’s function Gc{E) = Ac(E) — 
i'KUc{E) can be computed explicitly 



= E - Tl) ©(Tl - Mil). (4b) 

n 

where in Eq. (4b) the summation ji runs over the low (—) 
and high (+) energy bands of the unperturbed system 
and Ac = 0 for dimer, Ac = 1 for non-dimer sites. We 
can rewrite uc(E) as vc{E) = P^(E) pQ (E). Here, 

Qq (E) is the BLG unperturbed density of states per atom 
and spin for the jith energy band, 

^o(E) = Q(-P - Ul) Q(TI - Mil), (5) 

and the “projection factor” 


where M denotes the number of unit cells of the system, 
yields for each k four spin-degenerate eigenenergies [6, 7] 


Sa(k) = I (/lii + /i? + 4ig|/(k)|2 j , 
with the corresponding eigenstates 


( 2 ) 




The subscript a labels the conduction (a = +) and va¬ 
lence {a = —) bands, and its superscript counterpart fi 
distinguishes the high {/i = +) and low (// = —) en¬ 
ergy bands with respect to the charge neutrality point 
(at zero energy with Fermi momenta ±K). As in SLG, 
/(k) is the structural tight-binding function, /(k) = 
1 -h -h where ai and 3.2 are the primitive 

BLG lattice vectors, |ai| = |a 2 | = a = 2.46 A. 


SCATTERING FORMALISM FOR BILAYER 
GRAPHENE WITH A LOCAL IMPURITY 

To investigate scattering on a local perturbation V = 
E|c)(c| residing on the (non)dimer carbon site C, we em¬ 
ploy the projected retarded Green’s function Gc(E+) = 


specifies the relative contribution of the atomic orbital 
|c) to a group of /^-band eigenstates with energy E. The 
effective bandwidth D of the linearized theory is fixed by 
the consistency condition J^^[qq(E) + Qq (E)]dE = 1 
and equals VxTSrto — 6eV. 


BROADENING EFFECTS: REVERSING 
SPIN-RELAXATION TRENDS IN 
SINGLE-LAYER GRAPHENE 

The spin-relaxation in single-layer graphene (SLG), as 
a function of the carrier density, shows the opposite trend 
to that of bilayer graphene (BLG), despite the fact that 
the calculated bare 1 /r^ at zero temperature and no en¬ 
ergy broadening due to electron-hole puddles look simi¬ 
lar. As we argue in the manuscript, the measured l/r^ 
of SLG is strongly affected by both temperature and en¬ 
ergy broadening, reversing the bare trend. Temperature 
broadening is more important in SLG than in BLG be¬ 
cause the singlet-triplet resonance peaks in SLG are less 
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FIG. 1. Calculated spin-relaxation rate of ultraclean SLG 
(no energy broadening) at temperatures 0, 77, and 300 K. 
The magnetic adatom concentration is 1 ppm. While at 0 
and 77 K one sees the bare resonance trend, at 300 K the 
dependence of I/ts versus energy is reversed at low energies. 


separated than in BLG. To guide potential experiments, 
we show our calculations of l/r^ for SLG at different 
temperatures, and no energy broadening, in Fig. 1. Pa¬ 
rameters are the same as in Ref. [8]. While the bare 
trends, reminding of the resonance scattering, are visi¬ 
ble still at 77 K, at room temperature we see a complete 
reversal of the trend. This means, that at room temper¬ 
ature, even in ultraclean graphene with no electron-hole 
puddles, we do not expect to see the characteristic bare 
resonance shape. Instead, the experiment would see the 
usual “Elliot-Yafet” characteristics. 

However, at lower temperatures, say below 100 K, ex¬ 
periments could probe the spin-relaxation dependence on 
the electron density, by making the samples cleaner. In 
Fig. 2 we plot the calculated spin-relaxation rate at 
77 K at different electron-hole puddle broadening lev¬ 
els. At low broadenings, below 20 meV, the resonance 
trend should be visible, and SLG and BLG should exhibit 
the same “Dyakonov-Perel” trends in the spin-relaxation 
rate, if additional broadening due to different adatom 
types would be absent. The broadening starts to reverse 
the trend at 30 meV, at which I/r^ is rather flat at low 


energies, completing the reversal above 40 meV of broad¬ 
ening. We checked that the same happens at 4 K—also 
here the trend reverses for 40 meV of broadening. 
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FIG. 2. Galculated spin-relaxation rate of SLG at 77 K. The 
magnetic adatom concentration is I ppm. At 10 and 20 meV 
broadening, the spin-relaxation rate I/r^ exhibits a resonance 
trend at low energies, reminiscent of the “Dyakonov-Perel” 
shape, as also observed in BLG. At 30 meV the rate becomes 
flat, while at 40 meV the reversal of the trend towards the 
“Elliott-Yafet” shape takes place. 
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